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ABSTRACT 

Measurements  of  forward  scattered  light  were  used  to 
determine  the  particulate  behavior  in  the  combustor  and 
across  the  exhaust  nozzle  of  a  small  solid  propellant  rocket 
motor.  A  MALVERN  2600  was  successful  in  measuring 
multimodal  particle  distributions,  and  a  locally  designed 
system  was  successfully  applied  to  monomodal  distributions. 
The  propellant  tested  was  a  reduced  smoke  propellant 
containing  one-percent  zirconium  carbide  by  weight.  The 
results  suggested  that  ZrC  reacts  to  completion  on  or  very 
near  the  propellant  surface.  Zirconium  oxide  (Zr02) 
agglomerates  are  believed  to  form  on  the  surface.  Surface 
agglomerate  size  decreased  with  increasing  pressures  to 
approximately  500  psig.  Agglomerates  broke  up  in  the  nozzle 
convergence  to  approximately  18  microns,  independent  of  the 
pressure  or  nozzle  entrance  particle  size.-  Particle 
breakup  occurred  in  the  exhaust  plume..  System  limitations, 
sources  of  error  and  suggested  improvements  are  discussed. 
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I.   INTRODUCTION 

The  use  of  metallized  propellants  offers  the  advantage 
of  improved  solid  rocket  motor  performance  in  terms  of 
increased  energy,  specific  impulse,  and  suppression  of  high 
frequency  acoustic  combustion  instabilities.  Aluminum  is 
currently  the  primary  metallic  additive  of  choice  when 
exhaust  plume  visibility  is  not  a  mission  consideration. 
However,  the  requirement  for  reduced  smoke  propellants  for 
use  in  tactical  air-to-air  applications  has  dictated  the  use 
of  alternate  high  energy  propellants  which  utilize  only 
small  amounts  of  additives  such  as  zirconium,  zirconium 
carbide,  and  aluminum  oxide  to  provide  acoustic 
stabilization. 

The  improvements  achieved  via  such  additives  are  not 
without  penalty,  certain  loss  mechanisms  are  introduced. 
The  specific  impulse  efficiency  of  a  propellant  containing 
metal  is  generally  less  than  that  of  a  base  propellant  due 
to  incomplete  combustion  of  the  metal  in  the  motor  and  the 
presence  of  condensed  metal  and  metal  oxides  in  the  nozzle 
flow.  Particle-gas  velocity  lag  represents  the  primary 
source  of  performance  loss  within  the  exhaust  nozzle. 
Thermal  lag  generally  has  less  effect  on  performance.  Yet, 
neglecting  thermal  and  velocity  lag  as  loss  mechanisms,  the 
simple  fact  that  the  condensed  particulates  do  not  expand  as 


the  other  combustion  gases  expand,  results  in  a  reduction  of 
available  specific  impulse. 

Netzer  et  al.  [Ref.  1],  point  out  that  the  behavior  of 
metal  as  it  leaves  the  burning  surface  is  not  a  well-defined 
phenomenon.  Some  particles  upon  reaching  the  burning 
surface  depart  immediately,  while  others  agglomerate  on  the 
surface  before  passing  into  the  gas  flow.  The  mechanics  of 
metal  droplet  combustion  as  discussed  by  Schorr  [Ref.  2], 
indicate  that  resultant  particles  in  the  range  of  1-2 
microns  and  smaller  are  of  significance  in  determining 
exhaust  plume  signature,  while  those  particulates  of  5 
microns  and  greater  are  important  in  the  determination  of 
two-phase  flow  losses.  The  problem  of  predicting 
performance  depends  on  an  accurate  determination  of  these 
two-phase  flow  losses  which  are  a  strong  function  of 
particle  size.  The  aluminum  oxide  that  results  from  the  gas 
phase  combustion  of  aluminum  is  generally  less  than  two 
microns  in  size.  These  small  particles  have  little  velocity 
lag.   Thus,  making  them  smaller  would  serve  little  benefit. 

Several  complex  computer  codes  attempt  to  model  the 
combustion  process  in  order  to  predict  solid  rocket  motor 
performance.  The  SPP  (Solid  Performance  Program)  [Refs. 
3,4],  is  one  such  code.  These  programs  are  semi-empirical 
in  nature  and  require  an  a  priori  knowledge  of  particle  size 
distribution  entering  the  exhaust  nozzle.   This  distribution 


must  be  experimentally  determined,  as  the  problem  is 
currently  intractable  with  regard  to  analytic  methods. 

The  following  experimental  methods  have  been  employed  in 
previous  work  at  the  Naval  Postgraduate  School: 

(a)  collection  of  exhaust  products 

(b)  light  transmission  measurements 

(c)  light  scattering  measurements 

(d)  holography  and  high  speed  motion  pictures. 

A  summary  of  previous  research  conducted  at  the  Naval 
Postgraduate  School  combustion  laboratory  is  contained  in 
Reference  1.  The  present  investigation  focused  on  the  use 
of  light  scattering  techniques  to  determine  particle  size. 
The  method  is  non-intrusive  and  suited  to  measurements  in 
both  the  motor  and  exhaust  nozzle  environments.  The  speed 
of  the  flow  precludes  successful  application  of  most 
holographic  techniques  within  the  exhaust  nozzle. 

Heretofore  most  studies  have  been  directed  at  aluminized 
composites,  and  much  is  known  about  the  combustion  of 
aluminum.  The  behavior  of  other  metallic  additives  have  not 
been  as  fully  explored.  For  example,  zirconium  carbide  is 
used  as  a  stability  additive.  However,  it  is  not  known 
whether  surface  agglomeration  occurs,  or  whether  the 
zirconium  burns  to  zirconium  oxide.  Both  of  these  processes 
will  affect  the  size  distribution  and  density  of  the 
particles  in  the  gas  phase.    These  in  turn  affect  the 


ability  of  the  particles  to  provide  damping  for  a  specific 
frequency  of  pressure  oscillation. 

The  goals  of  this  investigation  were: 

(1)  to  refine  previously  developed  experimental  and 
diagnostic  techniques,  and 

(2)  to  use  the  apparatus  to  determine  the  effects 
operating  pressure  and  position  within  the  motor  and 
exhaust  nozzle  on  the  particle  size  distributions 
from  a  reduced  smoke  propel lant. 

The  propellant  studied  was  a  low  smoke,  composite  which 

contained  one  percent  zirconium  carbide  by  weight.    The 

behavior  of  particulates  at  various  positions  in  the  motor 

and  exhaust  nozzle  was  examined  as  a  function  of  chamber 

pressure. 


II.   THEORETICAL  BACKGROUND 

The  performance  of  a  solid  rocket  motor  utilizing 
metallized  propellants  is  significantly  dependent  on  the 
diameter  of  the  condensed  particulates  present  in  the  flow. 
A  non-intrusive,  in-situ  method  is  the  most  desirable  means 
to  determine  particle  size,  since  the  measurement  is  real 
time  and  the  environmental  conditions  (pressure, 
temperature,  gas  composition,  etc.)  are  undisturbed.  The 
method  chosen  for  the  present  research  was  the  dif fractively 
scattered  light  technique.  Very  large  particles  cannot  be 
measured  with  this  technique  due  to  the  transmitted  light 
interference  with  the  light  which  is  scattered  at  very  small 
forward  angles.  Very  small  particles  scatter  at  larger 
angles,  but  the  intensity  is  low,  making  detection 
difficult.  Absorbtion  and  refraction  also  become  important 
for  very  small  particles.  Particle  concentrations  are  also 
limited  to  those  which  do  not  exhibit  significant  multiple 
scattering  effects.  Within  these  limits,  the  mean  particle 
size  is  readily  determined  from  diffraction  techniques,  but 
size  distributions  cannot  easily  be  determined.  In  general, 
a  distribution  type  must  be  selected  a  priori,  i.e.,  Rosin- 
Rammler,  log-normal  or  other  such  distribution  must  be 
assumed.  Yet,  the  diffraction  technique  is  non-intrusive 
and   can   be   used   in   the   internal   motor   and   nozzle 


environment,  where  velocities  and  particle  concentrations 
can  be  high. 

The  scattering  and  absorbtion  of  a  plane  monochromatic 
electromagnetic  wave  by  an  isotropic  homogeneous  sphere  is 
described  by  Mie's  solution  to  Maxwell's  equations  [Ref.  5], 
and  the  general  theory  of  scattering  developed  by  Mie  is 
presented  by  Van  de  Hulst  [Ref.  6].  According  to  Mie's 
theory,  the  characteristics  of  scattered  light  depend  upon 
the  particle  size,  shape,  internal  structure,  refractive 
index  with  respect  to  the  surrounding  medium,  the  scattering 
angle  and  the  incident  wave  polarization  and  wavelength 
[Ref.  7].  The  solution  to  Maxwell's  equations  for  multiple 
scattering  in  a  concentrated  sample,  where  interaction  of 
scattering  centers  is  significant,  presently  represents  an 
intractable  problem,  even  if  boundary  conditions  could  be 
specified.  Even  the  solution  without  multiple  scattering  is 
quite  complex  and  difficult  to  incorporate  into  diagnostic 
equipment. 

When  the  particles  are  much  larger  or  much  smaller  than 
the  wavelength  of  the  incident  light,  the  equations  simplify 
to  Fraunhofer  diffraction  and  Rayleigh  scattering, 
respectively.  The  light  scattered  from  particles  within 
these  restricted  ranges  can  be  measured  and  used  to  readily 
determine  particle  size.  The  size  of  particles  of  interest 
in  solid  propellant  combustion  depend  on  the  application. 
When  considering  two  phase  flow  losses,  particle  sizes  of 


interest  are  generally  five  microns  or  greater.  Jones  [Ref. 
8],  has  shown  that  Fraunhofer  diffraction  can  be  used  for 
particles  as  small  as  one  micron  when  a  He-Ne  laser  is  used. 
Fraunhofer  diffraction  is  applicable  in  the  case  of  larger 
particles  because  the  forward  scattered  light  in  the  center 
lobe  is  due  principally  to  diffraction.  The  mechanisms  of 
reflection  and  absorbtion  do  not  significantly  contribute  to 
the  center  lobe  in  this  instance  (this  was  shown  by 
Hodkinson  [Ref.  9]).  If  these  latter  effects  are  present  to 
a  significant  extent,  then  the  complete  Mie  equations  must 
be  used. 

The  classical  results  of  Fraunhofer  diffraction  from 
particles,  i.e.,  the  discrete  rings  associated  with  a 
monodispersion,  are  typically  not  observed  in  practical 
combustion  experiments.  Typically,  the  particulates  are 
polydispersed  and,  in  many  cases,  a  multimodal 
polydispersion  is  present.  The  intensity  distribution  due 
to  a  monodispersion  is  given  by  Fraunhofer  diffraction  [Ref. 
10]  as, 


2J  (a  sin  6)  2 

*i§i   =   [— i . s-]  (D 

1(0)  a  sin  9 


where  J]_  is  the  first  order  Bessel  function,  a  =  tt  D/X  ,  9  is 
the  scattering  angle  for  the  forward  direction,  D  is  the 
particle  diameter,  and  A  is  the  wavelength  of  the  incident 


light.     When  considering   small   angles   (9  <   15°),   the 
expression  simplifies  to, 
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Dobbins  [Ref.  11]  showed  that  the  normalized  integrated 
intensity  of  forward  scattered  light  from  a  polydispersion 
of  large  particles  is  given  as, 
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where  Nr(D)  is  a  distribution  function  describing  the 
proportion  of  particles  with  diameter  D  in  the  sample. 
Dobbins  et  al.  [Ref.  11]  used  the  Upper  Limit  Distribution 
Function  (ULDF)  developed  by  Mugele  and  Evans  [Ref.  12]  and 
the  same  approach  was  followed  in  a  portion  of  this 
investigation.  The  volume-to-surface  mean  diameter  (Sauter 
mean  diameter)  D32  of  a  polydispersion  is  defined  by, 
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The  significance  of  D32  (versus  D43,  etc.)  is  that  the 
illumination  profiles  for  different  size  distributions  of 
the  ULDF  plotted  against  (tD326//)  are  brought  into  near 
coincidence  [Ref.  11]. 

Buchele  [Ref.  10]  showed  that  for  TTD326/A  less  than 
three,  a  gaussian  curve  can  be  used  for  the  scattered  light 
intensity  which  closely  matches  the  theoretical  intensity 
profile  for  a  polydispersion  obtained  by  integrating 
Equation  (3).  The  gaussian  expression  presented  by  Buchele 
is  given  as, 


m|  -  exp[-(.57^D32e/X)2]  (5) 


The  data  reduction  technique  developed  by  Netzer  and  Harris 
[Refs.  1/13]  is  based  on  Buchele' s  work.  The  method 
involves  comparison  of  the  experimentally  obtained  intensity 
profile  with  that  produced  by  Equation  (5) .  Use  of  Equation 
(5)  requires  the  value  of  1(0),  the  forward  scattered 
intensity  at  theta  equal  zero.  Since  transmitted  light 
dominates  at  very  low  scattering  angles,  1(0)  cannot  be 
measured.  A  ratioing  technique  eliminates  this  problem. 
Equation  (5)  can  be  used  to  ratio  I(e2)/I(6i)  as  follows, 

TTfJ    =  exP{-D3V^)2(e2  -  el)}  .         (6) 


This  eliminates  the  need  to  determine  1(0).  However,  I(5]_) 
at  the  "minimum"  forward  angle  9]_  must  be  taken  from  the 
experimentally  measured  profile  at  a  location  where  the 
transmitted  light  effects  are  not  significantly  present. 
The  resultant  profile  I(e2)/I(9]_)  versus  ©2  can  then  be 
generated  for  a  specified  D32.  D32  is  varied  interactively 
until  the  profile  generated  by  Equation  (6)  matches  the 
measured  profile.  This  procedure  yields  D32.  However,  this 
approach  has  an  obvious  disadvantage  in  that  it  is  only 
accurate  if  the  polydispersion  is  monomodal. 

A  Malvern  2600  particle  sizer  was  also  used  in  the 
present  investigation.  The  Malvern  software  is  similarly 
based  upon  Fraunhofer  diffraction.  However,  it  is  based  on 
the  inversion  of  the  measured  scattered  data  to  obtain  a 
particle  size  distribution.  The  Malvern  allows  use  of  a 
selected  type  of  particle  size  distribution  (Rosin  Rammler, 
log-normal,  normal)  or  of  a  model  independent  mode.  The 
model  independent  method  is  described  by  the  Malvern  2  600 
Reference  Manual  [Ref.  14]  as  a  constrained  non-linear  least 
squares  method,  which  is  able  to  measure  multi-modal 
particle  distributions  with  high  resolution.  The  obvious 
advantages  of  the  Malvern  2600  are  the  ability  to  measure 
multi-modal  particle  size  distributions,  and  it  has  been 
independently  tested.  Dodge  et  al.  [Ref.  15]  conducted  a 
comparison  of  the  Malvern  device  and  a  phase  doppler  system 
and   found  that  there  was  excellent  agreement  in  point 
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measured  Sauter  mean  diameters  between  the  deconvoluted 
laser  diffraction  data  and  the  phase  doppler  data.  The 
Malvern  system  is  currently  one  of  the  most  widely  used 
laser  diffraction  devices  and  it  shows  great  promise  for 
application  to  particle  sizing  in  solid  rocket  motors. 

Other  work  has  been  done  in  the  area  of  analytic  methods 
that  appears  to  provide  the  basis  for  refinements  in 
particle  sizing  and  determination  of  particle  size 
distributions.  Koo  and  Hirleman  [Ref.  16]  evaluated  the  use 
of  integral  transform  technigues  on  Fraunhofer  diffraction 
patterns.  These  technigues  offer  the  advantage  of  direct 
inversion  of  scattering  data  to  obtain  a  size  distribution. 
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III.       EXPERIMENTAL   APPARATUS 

A.  INTRODUCTION 

Laser  diffraction  systems  have  become  widely  used  in  the 
field  of  particle  and  droplet  size  measurements.  Fraunhofer 
diffraction  is  the  underlying  principle  involved.  Hodkinson 
[Ref.  9]  demonstrated  that  for  particles  greater  than 
several  microns  in  size,  the  primary  contribution  in  the 
near  forward  small  angle  region  is  from  diffraction  rather 
than  reflection  or  refraction  and  Jones  [Ref.  8]  has  proved 
that  in  this  same  region,  Fraunhofer  diffraction  may  be 
applied  to  particles  as  small  as  one  micron  when  a  helium- 
neon  laser  is  chosen  as  the  source  of  illumination.  Figure 
3.1  details  a  simple  schematic  for  a  laser  light  diffraction 
particle  sizer.  Such  a  system  typically  consists  of  a  He-Ne 
laser  source,  beam  expander  (spatial  filter  and  collimating 
lens) ,  Fourier  transform  lens,  and  a  detector. 

Typically,  small  particles  diffract  light  at  large 
angles  and  large  particles  diffract  light  at  small  angles. 
Figure  3.2  depicts  the  diffraction  process.  Note  that  the 
undiffracted  light  is  focused  to  a  single  point. 

B.  ROCKET  MOTOR 

The  series  of  experiments  conducted  required  the  use  of 
two  small  solid  propellant  rocket  motors.  The  motor 
depicted  in  Figure  3.3  is  the  same  as  that   used   by   Pruitt 
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La/ a®  particles 
scatter  at  low  angles 


Detector  measures 
Integral  scattering  of  all 
particles  simultaneously. 


Central 
A  Detector 


Small  panicles 
scatter  at  high  angles 


Radii 


Figure  3.2   Scattering  Properties  of  Particles  as 
a  Function  of  Size  [P.ef.  14] 
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Figure  3 . 3   Motor  Components 


Figure  3.4   Short  Motor  Components 


15 


[Ref.  17].  The  chamber  diameter  was  2.00  inches  and  the 
propellant  was  cylindrical,  1.98  inches  in  diameter  and 
1 .  00-1.  25. -inches  thick.  All  tests  conducted  were  done  with 
end  burning  grains.  The  long  motor  (Fig.  3.3)  was  used  to 
obtain  in-situ  measurements  of  particle  size  prior  to  the 
nozzle  entrance  and  in  the  exhaust.  Figure  3.4  depicts  the 
short  motor  which  was  used  to  obtain  in-situ  particle  size 
measurements  just  aft  of  the  propellant  burning  surface. 
This  particular  motor  had  been  used  previously  in  connection 
with  holographic  imaging.  For  additional  specifications  see 
Reference  18. 

Both  motors  were  windowed,  and  the  windows  were  kept 
clean  using  a  nitrogen  purge.  The  nitrogen  was  "spread" 
using  a  sintered  metal  sleeve  which  was  inserted  into  the 
aperture. 

Ignition  was  accomplished  using  a  BKNO3  ignitor  which 
was  fired  by  means  of  a  nichrome  filament  energized  by  a  12 
VDC  source. 

The  propellant  was  bonded  into  the  motor  casing  with  a 
self-vulcanizing  silicone  rubber  (RTV) .  The  propellant 
studied  was  a  reduced  smoke  propellant  with  87%  solids 
loading.  Table  3.1  presents  the  major  propellant  ingredient 
weight  percentages. 

C.   LASER  LIGHT  SCATTERING  APPARATUS 

The  locally  designed  laser  diffraction  system  makes  use 
of  an  8  mWatt  He-Ne  source  in  the  same  general  configuration 
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TABLE  3 . 1 
PROPELLANT  INGREDIENTS 

Ingredient  Approx.  Wt.  % 
R4  5M  9.2 

RDX  4  .  0 

AP  8  2.0 

ZrC  1.0 

depicted  in  Figure  3.1.  (A  complete  description  of  the 
apparatus  is  contained  in  Reference  2  0  [see  Figures.  3.5  and 
3.6]).  The  scattered  light  is  passed  through  a  wide-pass 
laser  line  filter  and  a  500  mm  focal  length  Fourier  trans- 
form lens,  which  focuses  the  light  onto  the  detector.  The 
detector  is  a  linear  photodiode  array  (Reticon  RL  1024G) 
which  has  1024  photodiodes  on  25  micron  centers.  The  array 
is  mounted  on  a  "mother"  board  which  contains  the  sample- 
and-hold  circuits.  Each  photodiode  measures  the  intensity 
of  the  scattered  light  and  outputs  a  proportional  voltage. 
The  array  is  scanned  at  a  rate  of  3  3  kHz  and  a  6  msec  blank- 
ing pulse  is  used  to  sequence  the  sweeps.  The  range  of 
scattering  angles  within  the  field  of  view  is  determined  by 
a  number  of  factors.  The  upper  limit  of  scattering  angle  is 
determined  by  the  diameter  and  focal  length  of  the  focusing 
lens,  the  distance  between  the  focusing  lens  and  the 
particles,  the  diameter  of  the  motor  window,  and  the 
height/position  of  the  diode  array  [Ref.  1] .   The  range  of 
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scattering  angles  for  this  investigation  was  .01-. 06 
radians,  corresponding  to  an  offset  of  5  mm  of  the  array  and 
.02-. 07  radians  corresponding  to  a  10  mm  offset.  Because  of 
the  very  small  detectors,  the  system  is  not  capable  of 
measuring  very  weak  scattering  intensities. 

The  second  system  utilized  was  a  commercially  available 
particle  size  analyzer  manufactured  by  MALVERN  Instruments. 
The  particular  instrument  used  in  this  series  of  experiments 
was  the  MALVERN  2600  HSD.  Figure  3.7  depicts  the  MALVERN 
2600  system,  which  utilizes  a  2  mW  He-Ne  laser  with  a 
collimated  beam  diameter  of  9  mm  (Table  3.2  details  laser 
specifications) .  One  of  the  principle  advantages  of  the 
MALVERN  system  is  its  use  of  an  annular  photodiode  array. 
The  array  is  a  monolithic  detector  which  consists  of  31 
rings.  It  has  increasing  thicknesses  of  the  annular  detec- 
tor elements  with  increasing  circumference,  resulting  in  an 
increase  of  detector  area  as  radius  increases.  Ring  arrays 
permit  much  weaker  intensities  to  be  measured  than  can  be 
accomplished  with  linear  arrays.  The  dynamic  range  of  the 
instrument  is  180:1  on  any  of  the  three  range  settings,  with 
an  advertized  accuracy  of  plus  or  minus  4%  of  volume  mean 
diameter  [Ref.  19].  The  detector  may  be  used  in  conjunction 
with  any  one  of  three  interchangeable  Fourier  transform 
lens,  with  focal  lengths  of  63,  100  and  300  mm  [Ref.  14]. 
The  respective  size  classes  associated  with  each  lens  are 
1.2-118,  1.9-188   and   5.8-564   microns,   with   a   sub-class 
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Figure  3.6   Hewlett  Packard  HP  9836S  Computer  and 
HP  6942A  Multiprogrammer 


20 


Receiver  Unit 


Transmitter  Unit 


Analyser  Beam 


Range  Lens 


Beam  Expander 


Optical  Bed 


2600  Optical   Unit   -  Components  &   Controls 


Olivetti  M24  AT&T  PC 


Olivetti  M21 


Receiver  Unit 


Oonputer  to  Optical  Unit  Connection 

Figure    3.7      MALVERN    2600    Particle    Sizer    [Ref.    14] 
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TABLE  3.2 
LASER  SPECIFICATIONS 

Helium-Neon  Laser  (Locally  Designed  System) 

1.  Manufacturer:  Spectra-Physics 

2.  Model:  147 

3.  Type:  He-ne  Class  111B 

4.  Output  Power:  8  mWatt 

5.  Beam  Diameter:  .92  mm 

6.  Beam  Divergence:  .87  mrad 

Helium-Neon  Laser  (MALVERN  System) 

1.  Manufacturer:  MALVERN  Instruments 

2.  Model:  2600 

3.  Type:  He-Ne  Class  111B 

4.  Output  Power:  2  mWatt 

5.  Beam  Diameter:  9  mm 

6.  Beam  Divergence:  Not  Available 


which  extends  down  to  .5  microns  [Ref.  14].   The  100  mm  lens 
was  used  for  this  set  of  experiments. 

D.   DATA  ACQUISITION  AND  REDUCTION 

The  data  acquisition  system  associated  with  the  locally 
designed  laser  diffraction  system  consisted  of  a  Hewlett 
Packard  6942A  multiprogrammer  which  was  controlled  by  an  HP 
9836S  computer.  The  data  acquisition  program  developed  by 
Harris  [Ref.  13]  was  the  system  driver  in  the  acquisition 
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mode.  In  the  course  of  a  motor  firing  sequence,  a 
background  measurement  was  recorded  prior  to  ignition.  This 
mean  intensity  profile  was  then  subtracted  from  the  mean 
intensity  profile  measured  with  particles  present.  This 
procedure  accounts  for  noise  present  due  to  ambient  light 
and  corrects  for  the  characteristics  of  individual  diodes. 
The  raw  data  were  displayed  on  the  HP  98  3  6S  CRT  so  that  any 
erroneous  scans  could  be  excluded  prior  to  data  reduction. 

A  symmetric  moving  average  type  digital  filter  was  then 
applied  to  the  profile  in  order  to  smooth  the  curve  [Ref. 
1]  .  This  type  filter  does  not  have  the  phase  lag  that  is 
inherent  in  analog  filters.  The  data  reduction  algorithm  is 
based  on  the  work  done  by  Buchelle  [Ref.  10],  and  is  fully 
described  in  Reference  20.  The  chief  disadvantage  of  the 
method  lies  in  the  fact  that  there  is  no  way  to  account  for 
the  actual  distribution  of  the  particulates.  The  Sauter 
mean  diameter  (SMD)  determined  by  this  algorithm  will  be 
that  which  is  associated  with  the  dominant  mode  (within  the 
field  of  view  of  the  array)  of  a  multimodal  polydispersed 
sample. 

The  MALVERN  system  represents  a  vast  improvement,  in 
that  not  only  is  SMD  determined,  but  the  volume  distribution 
is  characterized  as  well.  The  MALVERN  2600  data  acquisition 
system  is  centered  about  an  AT&T  63  00  PC.  The  system  driver 
is  the  MALVERN  Master-Sizer  operating  system.  This  version 
of  the  MALVERN  software  uses  a  command-based  language  that 
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allows  the  user  to  specify  the  pertinent  parameters 
associated  with  the  operation  of  the  instrument.  The 
Master-Sizer  program  allows  the  user  to  select  from  any  one 
of  four  models.  The  model  types  available  are  Rosin- 
Rammler,  log-normal,  normal,  and  model  independent.  The 
model  independent  option  of  the  Master-Sizer  program  is  a 
process  of  constrained  least  squares  fitting  of  theoretical 
scattering  characteristics  to  observed  data  [Ref.  14].  It 
is  usual  to  allow  the  best  fit  result  to  be  obtained  with  no 
assumed  form  of  size  distribution  and  this  allows  the 
characterization  of  polydispersed  multimodal  distributions 
with  high  resolution  [Ref.  14].  The  capability  of 
characterizing  multimodal  distributions  is  invaluable  in 
work  dealing  with  the  combustion  of  metallized  propellants, 
in  that  the  distribution  of  particulates  is  most  often 
multimodal.  It  is  for  this  reason  that  the  model 
independent  mode  was  used  exclusively  throughout  this 
investigation . 
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IV.   RESULTS  AND  DISCUSSION 

A.   INTRODUCTION 

This  investigation  was  based  on  three  major  objectives: 

(1)  refine  previously  developed  experimental  and 
diagnostic  techniques  by  identifying  specific  causal 
factors  of  measurement  error, 

(2)  incorporate  the  use  of  a  commercial  particle  size 
analyzer  (MALVERN  2600)  to  validate  and  quantify  the 
operating  limits  of  the  locally  designed  laser 
diffraction  system  and, 

(3)  use  the  apparatus  to  determine  the  effects  of 
operating  pressure  and  position  within  the  motor  and 
exhaust  nozzle  on  particle  size  distribution 
associated  with  the  combustion  of  a  reduced  smoke 
propellant . 

Previous  investigations  at  NPS  had  been  based  on  the  method 

of  Buchele  [Ref.  10]  which  yields  Sauter  mean  diameter,  but 

does   not   provide   any   information   as   to   the   actual 

distribution  of  particulate  size.   The  introduction  of  the 

MALVERN  system  allowed  for  a  characterization  of  volume 

distribution.    The  common  parameter  between  the  locally 

designed  system  and  the  MALVERN  was  Sauter  mean  diameter 

(D32).   This  parameter  provided  the  basis  for  comparison  of 

the  two  diffraction  systems.   Both  systems  were  used  in  the 

initial  evaluation  of  the  reduced  smoke  propellant.   Once 

the  limitations  of  the  locally  designed  instrument  were 

established,  the  MALVERN  was  used  more  extensively  to  gather 

data  in  both  the  motor  and  exhaust.   The  distinct  advantage 

of  being  able  to  characterize  multimodal  particle  size 
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distributions  made  the  MALVERN  the  better  suited  instrument 
for  the  through  motor  measurements.  In  the  exhaust,  the 
annular  detector  of  the  MALVERN  system  allowed  for  the 
measurement  of  much  weaker  intensities  as  opposed  to  the 
locally  designed  system  (which  is  also  referred  to  as  the 
NPS  system) . 

B.   SYSTEM  CALIBRATION 

Among  the  methods  available  for  calibration  of  laser 
diffraction  instruments,  the  more  commonly  used  include  the 
calibration  reticle  and  use  of  commercially  prepared  samples 
of  known  diameter.  The  latter  approach  was  used  in  this 
investigation.  The  types  of  samples  used  included 
polystyrene  spheres,  glass  beads  and  aluminum  oxide  powders. 
In  each  case,  the  sample  material  was  suspended  in  a 
solution  of  distilled  water. 

The  locally  designed  system  was  calibrated  using  the 
polystyrene  spheres,  with  diameters  of  5.1,  9.6  and  2  0.0 
microns,  respectively.  Table  4.1  contains  the  results  of 
these  calibration  tests.  Figures  4.1  and  4.2  are 
representative  samples  of  calibration  results. 

The  MALVERN  2600  was  calibrated  using  three  different 
materials,  polystyrene  spheres,  glass  beads  and  an  aluminum 
oxide  powder.  Results  are  presented  in  Table  4.1. 
Calibrations  from  the  two  systems  were  in  good  agreement. 
System  alignment  was  completed  with  the  container  of  water 
in  the  optical  path  and  a  background  measurement  was   taken 
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prior  to  introducing  particles.  A  glass  rod  was  used  to 
stir  the  mixture  to  ensure  that  the  particulates  remained  in 
suspension  during  the  measurement  process. 

A  primary  concern  early  in  the  investigation  was  to 
validate  the  ability  of  the  MALVERN  to  accurately 
characterize  a  multimodal  distribution  when  the  system  was 
operated  in  the  Model  Independent  mode.  In  order  to  test 
this  facet  of  the  Master-Sizer  program,  a  polydispersed 
bimodal  sample  consisting  of  5.1  and  20.0  micron  particles 
was  prepared.  The  results  are  depicted  as  a  volume 
distribution  in  Figure  4.3.  The  respective  modes  are 
centered  at  approximately  5.5  and  21.0  microns. 

When  multiple  scattering  becomes  significant,  the  mean 
particle  diameters  determined  from  single  scattering  theory 
are  too  small.   Gulder  [Ref.  281  has  determined  empirical 
corrections  for  the  Malvern  when  obscuration  is  greater  than 
0.5.   This  correction  was  utilized,  when  required,  in  the 
oresent  investigation. 

C.   RESULTS  FROM  REDUCED  SMOKE  PROPELLANT  EVALUATION 
1.   Test  Conditions 

The  series  of  experiments  conducted  covered  a  range 
of  pressures  from  approximately  70  to  450  psig.  The 
experiments  were  divided  into  three  pressure  classes,  150, 
250   and   400   psig.   At   each   class  of   chamber  pressure, 
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scattering  measurements  were  made  at  the  following 
locations: 

(a)  near  grain  (end  burning) 

(b)  prior  to  the  nozzle  entrance 

(c)  exhaust  (Pe  =  Pcritical) 

(d)  exhaust  (Pe  =  P0) . 

A  total  of  32  tests  were  completed,  24  of  which  produced 
usable  data.  Three  tests  failed  to  produce  data  because  the 
burst  disk  ruptured,  three  tests  of  the  short  motor  were 
invalidated  because  the  nitrogen  purge  system  failed  to  keep 
the  window  surfaces  clean,  one  test  failed  because  of  an  CD- 
ring  failure  in  the  motor  joint,  and  for  one  test,  data  were 
not  obtained  because  the  pressure  threshold  to  trigger  the 
acquisition  of  data  was  set  too  high. 

The  nozzles  were  designed/sized  for  isentropic  flow 
conditions  using  y  =  1.2.  The  ratio  of  specific  heats  was 
determined  using  PEPCODE,  an  equilibrium,  adiabatic 
combustion  code.  The  value  of  gamma  jotained  in  this  manner 
was  1.21.  Because  nozzle  efficiencies  should  be  quite  high 
(97-99%),  the  isentropic  assumption  was  made.  Deviations  in 
pressures  can  be  attributed  to  the  variations  in  ignition 
characteristics  and  propellant  properties  from  run  to  run, 
and  some  small  amount  of  nozzle  throat  erosion. 

The  propellant  composition  was  provided  by  the  Naval 
Weapons  Center  at  China  Lake,  California.  The  zirconium 
carbide  in  the  propellant  varied  in  size  from  2-20  microns, 
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with  a  D32  of  5.1  microns.   The  information  is  provided  here 
as   an   introduction   to   the   discussion   of   surface 
agglomeration  on  the  burning  grain. 
2 .   Zirconium  Carbide  Behavior 

The  size  of  agglomerates  that  form  on  the  surface  of 
burning  metallized  solid  propellants  are  determined  by  many 
factors,  including  metallic  powder  type,  propellant  type, 
size  of  original  metal  particulates,  concentration  of  metal, 
motor  chamber  pressure,  and  the  burning  rate  of  the 
propellant.  The  effect  of  chamber  pressure  on  agglomerate 
size  has  been  observed  to  be  guite  strong  for  aluminum.  As 
chamber  pressure  and  burning  rate  increase,  the  agglomerates 
become  smaller   [Refs.  21,22]. 

Zirconium  carbide  has  been  observed  to  burn 
vigorously  on  the  surface  of  reduced  smoke  propellants 
containing  AP  and  HTPB  [Ref.  23].  If  zirconium  carbide 
(ZrC)  burns  stoichiometrically  with  oxygen  it  results  in  a 
temperature  of  approximately  4400°K.  This  resultant 
temperature  is  much  greater  than  the  melting  temperature  of 
ZrC>2  (2988°K)  .  Table  4.2  delineates  the  temperatures 
relevant  to  the  discussion. 

A  plausible  behavior  of  the  ZrC  is  that  as  the 
material  is  exposed  at  the  grain  burning  surface,  oxygen 
(02)  from  the  ammonium  perchlorate  (AP)  flame  attacks  the 
carbon  to  produce  either  carbon  monoxide  (CO)  or  carbon 
dioxide   (C02) •   Oxygen   from  water  vapor   may  also   be 
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TABLE  4.2 
THERMOPHYSICAL  PROPERTIES  OF  Zr,  ZrC,  AND  Zr02 

Melting  Boiling 

Material  Temperature  (°K)       Temperature  (°K) 

Zr  2127  [Ref.  24]  4700  [Ref.  24] 

ZrC  3523  [Ref.  26] 

Zr02  2988  [Ref.  24]  4573  [Ref.  25] 

NOTE:   Adiabatic  flame  temperature  for  the  propellant  was 
determined  to  be  2916°K  and  the  adiabatic  flame 
temperature  for  the  reaction 

ZrC  +  202   "*"  Zr02  +  C02 

is  approximately  4400°K. 

involved.  The  oxygen  then  reacts  with  the  zirconium  to 
produce  a  porous  Zr02  surface  layer  [Ref.  24].  The  oxide 
should  be  a  liquid  since  the  local  temperature  at  the 
particle  surface  (as  a  result  of  the  surface  combustion  of 
ZrC  or  Zr  with  02)  should  be  much  greater  than  the  melting 
temperature  of  Zr02  (2988°K)  .  The  low  burning  rate 
associated  with  the  relatively  low  chamber  pressures  allows 
time  for  the  molten  Zr02  on  adjacent  particles  to  stick 
together  to  form  large  agglomerates  on  the  grain  burning 
surface,  prior  to  passing  into  the  flow.  At  higher  burning 
rates  (pressures)  the  particles  would  be  ejected  from  the 
surface  more  rapidly,  resulting  in  smaller  agglomerates. 
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3 .   Near-Propellant  Surface  Data 

The  test  results  obtained  near  the  grain  surface  are 
summarized  in  Table  4.3.  D32  values  are  plotted  in  Figure 
4.4.  As  can  be  seen  directly  from  these  results,  there  was 
a  strong  dependence  on  chamber  pressure  of  the  agglomerate 
size  produced  from  the  ZrC.  At  the  lower  pressures  (159  and 
188  psig) ,  the  distribution  of  particulates  ejected  from  the 
grain  surface  were  trimodal  and  bimodal,  respectively. 
Figures  4.5  and  4.6  detail  the  characterization  of  the 
particulate  distributions.  The  dominant  modes  were  centered 
at  approximately  30-37  microns  and  65  microns. 

The  Sauter  mean  diameters  determined  from  the 
scattering  data  were  6  3.5  microns  for  the  run  at  159  psig 
and  72.8  microns  for  the  run  at  188  psig.  This  was 
considered  to  be  good  repeatability,  run-to-run,  allowing 
for  variations  in  the  combustion  process  from  grain  to 
grain.  As  discussed  above,  many  of  the  ZrC  particles,  upon 
reaching  the  burning  surface,  apparently  react  to  form  Zr02 
and  agglomerate  before  being  ejected  into  the  gas  flow. 

Test  number  three  of  the  near  grain  results  (248 
psig)  showed  a  reduction  in  Sauter  mean  diameter,  with  a 
measured  D32  of  32.2  microns.  The  distribution  of 
particulates  at  this  pressure  was  still  essentially  bimodal 
as  depicted  by  Figure  4.7.  However,  the  dominant  modes  of 
the  distribution  were  now  centered  at  approximately  12-13 
microns  and  27  microns.   There  was  also  a   small   number   of 
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particulates  in  the  range  of  3-4  microns,  indicating  that 
some  of  the  particulates  at  this  pressure  (248  psig)  were 
ejected  directly  into  the  gas  stream.  The  448  psig  test 
resulted  in  a  D32  of  15.2  microns.  Figure  4.8  reveals  that 
the  distribution  of  particulates  continued  to  be  bimodal. 
It  is  apparent  from  Figure  4.8  that  a  much  larger  portion  of 
the  original  metal  (ZrC)  particles  was  ejected  directly  into 
the  gas  flow,  as  evidenced  by  the  well-defined  mode  centered 
between  3-4  microns.  The  general  observation  here  was  that 
as  motor  chamber  pressure  was  increased,  the  size  of 
agglomerates  on  the  grain  surface  decreased,  similar  to  the 
behavior  of  aluminized  propellants  reported  in  References  21 
and  22.  It  was  also  observed  that  a  greater  number  of  the 
original  ZrC  particulates  were  apparently  immediately 
ejected  into  the  gas  flow  as  they  reached  the  burning 
surface  as  chamber  pressure  was  increased.  In  all  cases, 
the  distributions  of  particulates  near  the  burning  surface 
were  bimodal  or  trimodal.  The  quantities  of  particles  in 
the  smallest  modes  of  the  trimodal  data  were  quite  small. 

The  most  difficult  obstacles  to  overcome  in 
obtaining  near  grain  light  scattering  measurements  were: 
(1)  keeping  the  window  surfaces  clean,  as  the  ignition 
process  (due  to  its  close  proximity  to  the  windows)  made  it 
difficult  to  effectively  purge  these  surfaces,  and  (2) 
multiple  scattering  effects  due  to  dense  concentrations  of 
particulate  matter  in  the  near  grain  region.   The  problem  of 
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multiple  scattering  has  been  addressed  above.   The  empirical 
corrections  presented  by  Gulder  [Ref.  28  ]  are  for  a  Malvern 
system  operated  in  the  model  independent  mode.   This  mode  was 
also  used  in  the  current  investigation.   When  obscuration  is 
approximately  0.99,  the  measured  D.,-  is  approximately  a  factor 
of  two  too  small.   Gomi  [Ref.  27  ]  has  also  discussed  the  effects 
of  multiple  scattering.   For  monodispersed  spheres  Netzer  and 
Powers  [Ref.  1  J  have  also  presented  data  for  correcting  D__ 
for  high  obscurations. 

The  more  difficult  problem  that  had  to  be  overcome 
in  the  near  grain  region  was  effectively  purging  the  window 
surfaces.  The  only  alternative  here  was  to  repeat  the 
experiment  until  usable  data  were  obtained.  As  noted 
earlier,  three  of  the  tests  conducted  in  the  near  grain 
region  were  invalidated  due  to  the  fact  that  the  window 
surfaces  had  become  contaminated  as  a  result  of  the  ignition 
process. 

4 .   Nozzle  Entrance  Data 

Table  4.4  presents  the  results  from  the  measurements 
obtained  prior  to  the  nozzle  entrance.  These  data  allow  one 
to  examine  the  change  in  particle  size  across  the  nozzle 
when  compared  to  the  exhaust  data  which  are  contained  in 
Tables  4.5  and  4.6.  The  D32  data  obtained  using  the  MALVERN 
are  also  plotted  in  Figure  4.16.  Figure  4.17  plots  these 
same  data  together  with  the  MALVERN  data  obtained   near   the 
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burning  surface.  Figures  4.9,  4.10  and  4.11  characterize 
the  volume  distributions  associated  with  the  150  psig  class 
of  experiments.  Note  in  Figures  4.9  and  4.10  that  the  dis- 
tributions remained  bimodal  in  nature,  with  the  dominant 
modes  being  centered  at  approximately  35  and  68-70  microns 
in  the  case  of  Figure  4.9  and  at  approximately  27  and  again 
at  68-70  microns  in  the  case  of  Figure  4.10.  In  the  case 
of  Figure  4.11,  note  that  the  profile  was  trimodal,  but 
shifted  toward  the  smaller  end,  with  dominant  modes  centered 
at  approximately  3-4  microns,  20  microns  and  at  38-40  mi- 
crons. The  Sauter  mean  diameter  obtained  from  the  scatter- 
ing data  was  30.7  microns  (Test  No.  3).  A  post-test 
alignment  check  done  with  the  MALVERN  system  in  the  align 
mode  revealed  that  some  significant  window  contamination  had 
occurred,  but  physical  examination  of  the  windows  showed 
very  little  evidence  of  contamination.  The  results  of  this 
one  test,  although  questionable,  are  included  for 
completeness.  Comparison  of  the  Sauter  mean  diameters  and 
volume  distribution  profiles  for  near  grain  and  nozzle 
entrance  measurements  completed  for  the  150  psig  class 
reveal  only  slight  differences.  This  occurred  even  though 
the  residence  time  (tres)  was  approximately  45  milliseconds. 
As  discussed  above  the  equilibrium  adiabatic  flame 
temperature  was  2916°K.  The  actual  temperature  may  have 
been  somewhat  less  than  this  due  primarily  to  heat  loss 
through  the  uninsulated,  stainless   steel  walls   of   the 
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combustor,  and  to  the  nitrogen  gas  used  to  purge  the  viewing 
windows.  The  observed  particle  behavior  would  be  expected 
if  the  ZrC  is  completely  consumed  near  (or  on)  the 
propellant  surface,  resulting  in  liquid  Zr02  particles 
flowing  in  a  low  velocity  (M  <<  .1)  environment. 

The  NPS  scattering  data  included  in  Table  4.4  will 
be  discussed  later  in  conjunction  with  the  comparison  of  the 
MALVERN  and  NPS  measurement  systems. 

Figures  4.12  and  4.13  represent  the  volume  distribu- 
tions associated  with  the  runs  completed  for  the  250  psig 
class.  This  represents  data  for  tests  8  and  9  of  Table  4.4. 
The  Sauter  mean  diameter  (SMD)  determined  from  scattering 
measurements  was  4  9.0  microns  corresponding  to  test  number  8 
and  48.0  microns  associated  with  test  number  9.  This  again 
represented  good  test-to-test  repeatability.  The  dominant 
modes  were  centered  at  approximately  the  same  locations  for 
both  tests.  Comparison  of  these  data  with  that  obtained  in 
the  near  grain  surface  region  showed  that  there  was  approxi- 
mately a  16  micron  particle  growth  over  the  length  of  the 
motor.  Note  that  only  one  measurement  was  taken  in  the  near 
grain  region.  If  the  point  was  valid,  then  particle 
collisions  and/or  wall  collisions/shedding  may  be  suspected. 

Test  numbers  10  and  11  represent  data  obtained  for 
the  400  psig  class  runs.  The  SMD's  determined  were  32.0  and 
27.6  microns,  respectively.  Figures  4.14  and  4.15  are  the 
volume   distribution   profiles   for  tests   10   and   11, 
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respectively.   Note  that  each  was  monomodal  in  nature.   By 
comparing  these  profiles  to  those  obtained  near  the  grain 
surface,  it  was  observed  that  the  3-4  micron  diameter 
particles  were  eliminated  as  they  passed  along  the  length 
of  the  combustor.   Again,  if  ZrC  is  assumed  to  react  to  Zr09 
on  or  near  the  propellant  surface,  particle  collisions  could 
explain  this  observation.   There  was  an  appreciable  increase 
in  particle  diameter  from  near  the  grain  surface  to  the 
nozzle  entrance  section.   Figure  4.16  is  a  graphical  depiction 
of  SMD  plotted  against  chamber  pressure.   Note  that  the  SMD 
decreases  with  increasing  chamber  pressure.   Figure  4.17 
presents  the  data  for  both  near-grain  and  nozzle-entrance 
measurements.   From  this  figure,  it  is  clear  that  particle 
growth/agglomeration  occurred  across  the  length  of  the 
combustor  for  pressures  above  200  psi. 
5 .   Exhaust  Nozzle  Data 

Table  4.5  is  a  listing  of  test  results  obtained  in 
the  exhaust  for  the  condition  of  nozzle  exit  plane  pressure 
equal  to  the  critical  pressure  at  the  nozzle  throat.   The 
nozzles  used  in  these  tests  had  a  45   converging  section  and 
a  short  flat  (.01"-. 015")  in  the  throat  section.   All  nozzles 
were  truncated  at  the  throat  after  the  .01"-. 015"  flat. 
Figures  4.18,  4.19,  4.20,  4.21  and  4.22  present  the  volume 
distributions  associated  with  tests  one  through  five,  respect- 
viely.    Figure   4.23   presents  the   results  obtained   on 
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test  6  using  the  NPS  system.  The  distance  downstream  of  the 
nozzle  exit  plane  varied  slightly  as  noted  in  Table  4.5  due 
to  the  geometry  of  the  apparatus  and  set-up.  All  the 
distributions  were  monomodal  with  the  exception  of  the  test 
conducted  at  the  lowest  pressure  (147  psig) .  Here  a  second 
broadly  distributed  mode  was  observed,  which  covered  a  range 
of  approximately  4.5-12.0  microns.  Except  for  this  point, 
the  SMD  varied  only  slightly  as  chamber  pressure  was  increased 
The  Weber  number  is  defined  as 

we  =  dag  Q?  (ug  "  uag) 2/ a 

where: 

dag  =  the  diameter  of  an  agglomerate  (m) 

Pq  =  the  density  of  the  gas 

uag  =  the  velocity  of  the  agglomerate 

Ug  =  the  velocity  of  the  gas 

g    =  the  surface  tension. 

It  is  the  ratio  of  inertial  forces  to  surface  tension 
forces.  Gany  [Ref.  22]  found  the  Weber  number  to  be  a  good 
correlation  parameter  in  the  analysis  of  the  break-up  of 
aluminum  oxide  agglomerates  in  the  nozzle  environment.  Gany 
[Ref.  22]  notes  that  for  low  Weber  numbers  the  droplets  are 
spherical,  and  tend  to  distort  when  the  Weber  number  exceeds 
about  four  (for  aluminum  oxide) .   The  distortion  continues 


64 


to  increase  as  Weber  number  increases  until  break-up  occurs 
for  Weber  numbers  in  the  range  of  12-20  [Ref.  22]. 

Figure  4.24  is  a  plot  of  SMD  versus  motor  chamber 
pressure.  It  is  apparent  from  the  curve  that  the 
agglomerates  broke  up  into  nearly  the  same  size  in  the 
converging  section  of  the  nozzle,  irrespective  of  chamber 
pressure  (unless  very  low)  or  particle  size  entering  the  nozzle 

The  equilibrium  adiabatic  combustion  temperature  was 
approximately  the  same  as  the  melting  temperature  for  Zr02. 
Thus,  even  though  some  thermal  lag  can  be  expected,  the 
particles  could  be  expected  to  rapidly  solidify  as  the  gas 
expansion  in  the  nozzle  occurs.  This  would  prevent  any 
further  reduction  in  particle  size. 

Table  4.6  contains  the  results  obtained  in  the 
exhaust  for  the  condition  of  nozzle  exit  pressure  equal  to 
ambient  pressure.  The  nozzles  all  had  converging  half 
angles  of  45°  and  diverging  half  angles  of  15°.  Nozzle 
throat  diameter  and  exit  diameter  were  determined  using 
isentropic  relations  for  a  ratio  of  specific  heats  equal  to 
1.2.  Figures  4.25  to  4.32  present  the  volume  distributions 
associated  with  tests  1  through  8,  respectively.  It  is 
interesting  to  note  that  all  the  distributions  were 
multimodal  with  the  exception  of  test  number  5,  which  was 
run  at  215  psig.  The  primary  factor  which  differentiated 
test  five  from  the  other  tests  was  the  distance  behind  the 
nozzle  exit  plane  at  which  the  measurements  were   obtained. 
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The  measurements  were  all  made  2.7  inches  downstream  of  the 
nozzle   exhaust   plane   except   for   test   five,   where   the 
distance  was  0.6  inches.   The  tests  run  with  the  2.7  inch 
distance  were  necessary  because  laser  diffraction  data  were 
being  obtained  simultaneously  through  the  motor,  and  the 
physical  dimensions  of  the  apparatus  limited  the  minimum 
distance  to  2.7  inches.   Test  5  results  indicated  that  the 
majority  of  the  particles  remained  monomodal  and  essentially 
unchanged  in  mean  diameter  (19  microns)  from  the  throat  to 
the  exit  plane.   The  measurements  made  further  aft  had  one 
mode  centered  at  the  same  19  micron  diameter.   However,  the 
number  of  these  particles  were  less,  and  significant  numbers 
of  smaller  particles  were  observed  in  a  bi-  or  tri-modal 
distribution.  The  reduced  particle  size  could  be  due  to  particle 
breakup.   However,  it  may  have  resulted  from  the  reduced 
number  of  particles  in  the  scattering  volume,  which  was  locat- 
ed through  the  center  of  the  spreading  exhaust  jet. 

The  presence  of  strong  thermal  gradients  along  the 
boundary  of  the  exhaust  jet  initially  caused  some  problem  in 
obtaining  exhaust  measurements.  During  initial  runs, 
scattering  data  indicated  the  presence  of  abnormally  large 
particulates  in  the  exhaust,  on  the  order  of  127  microns. 
Knowing  that  particle  size  entering  the  nozzle  was 
substantially  less  than  127  microns,  typically  around  70 
microns  for  150  psig  runs,  a  search  for  causal  factors  was 
pursued.    Conversation  with  Hughes  and  Taylor  at  MALVERN 
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Instruments  confirmed  the  idea  that  thermal  gradients  or 
density  gradients  were  indeed  biasing  the  scattered  light 
data,  i.e.,  the  presence  of  strong  thermal  gradients  along 
the  boundary  of  the  hot  gas  jet,  or  density  gradients 
associated  with  a  shock.  The  solution  to  this  difficulty 
was  a  rather  simple  procedure.  When  unstable  refractive 
index  of  the  suspending  medium  is  present  as  a  result  of 
strong  thermal  gradients,  a  steering  effect  occurs  that 
steers  light  to  low  angles,  hence  biasing  particle  size  to 
the  large  end  of  the  scale.  The  procedure  employed  when 
using  the  Master  Sizer  program  was  to  inhibit  data  from  the 
affected  inner  rings,  using  the  "KILLDATA"  command.  Once 
the  user  had  examined  the  raw  data  and  observed  which  rings 
had  been  affected,  then  those  rings  were  inhibited  by  the 
software  command  "KILLDATA."  Care  must  be  exercised  in  this 
procedure  to  ensure  that  no  scattering  data  will  be  lost  as 
a  result  of  inhibiting  data  from  the  affected  inner  rings. 
In  the  latter  case  there  is  no  recourse,  scattering  data  at 
that  point  becomes  no  longer  separable  from  the  steered 
light. 

Figure  4.33  is  a  graphical  depiction  of  particle 
size  (SMD)  versus  motor  chamber  pressure.   There  was   a 
slight  reduction  in  particle  size  with  increasing  pressure. 
Figure  4.34  plots  both  exhaust  conditions.  The  reduced  size 
in  the  downstream,  spreading  plume  is  apparent. 
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Exhaust  measurements  were  made  for  both  Pe=  p   and  Pe= 
pcritical  to  determine  if  any  change  in  particle  size  was 
incurred  as  a  result  of  the  supersonic  expansion  of  the  jet. 
Particle  breakup  downstream  of  the  nozzle  exhaust  was  observed. 
The  limited  data  at  approximately  200  psi  indicated  that  little 
breakup  occurred  within  the  nozzle,  downstream  of  the  throat. 
In  aluminized  propellants,  one  might  reasonably  expect  that 
vapor  condensation  occurring  during  the  expansion  process  might 
well  affect  the  size  distribution  of  exhaust  particles.   How- 
ever, in  the  case  of  zirconium  carbide,  the  temperature  in  the 
chamber  is  not  sufficiently  high  for  appreciable  amounts  of 
either  Zr,  ZrC  or  Zr02  to  be  present  in  a  vapor  phase. 
6 .   SEM  Analysis 

As  a  means  of  independently  verifying  the  scattering 
results,  residue  was  collected  from  the  face  of  the  nozzle 
converging  section.   This  residue  was  then  placed  in  a  solution 
of  lab  grade  acetone  and  put  into  an  ultrasonic  cleaner,  mixed, 
and  allowed  to  settle.   After  repeated  cleanings,  the  acetone 
was  evaporated  off  in  an  oven  under  reduced  pressure.   Using 
this  powder,  seven  samples  were  prepared  for  examination  with 
the  scanning  electron  microscope.   Three  of  the  samples  were 
prepared  using  a  non-conducting  epoxy  base  which  had  to  have 
a  gold  flashing  applied  in  order  to  provide  for  a  conducting 
surface.   The  remaining  samples  were  mounted  using  a  conduct- 
ing colloidal  carbon.   Figures  4.35  and  4.36  are  SEM  photo- 
graphs taken  of  the  samples.   Particles  were  counted  from  each 
sample  in  an  attempt  to  get  a  representative  size  distribution.  Some  218 
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5KU    MD  =  431111        S  =  00088    P  =  88883 


Figure  4.35   SEM  Photograph  of  Residue  Collected 
from  Exhaust  Nozzle 


Figure  4.36   SEM  Photograph  of  Residue  Collected 
from  Exhaust  Nozzle 
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particles  were  counted  with  sizes  ranging  from  1-20  microns. 
The  smaller  particles  (3  microns  and  less)  were  difficult  to 
see,  but  did  exist  in  the  sample.  No  particle  larger  than 
20  microns  was  observed,  but  this  by  no  means  implies  that 
no  larger  particles  were  present.  Collection  of  residues 
and/or  exhaust  products  is  a  crude  method.  When  collected 
in  small  quantities  from  wall  surfaces  it  cannot  yield  a 
valid  statistical  sample.  However,  it  does  provide  some 
reassuring  insight  as  to  what  was  physically  present  in  the 
sample. 

The  particles  were  grouped  strongly  in  2  bands,  one 
band  from  5-12  microns  and  the  other  from  14-17  microns. 
The  sample  was  taken  from  the  427  psig  run  where  the  mean 
particle  size  (SMD)  prior  to  the  nozzle  entrance  was  27.6 
microns  and  in  the  nozzle  throat  region  it  was  19 
microns.  Thus,  the  SEM  data  were  in  reasonably  good 
agreement  with  the  optical  data.  In  Figures  4.3  5  and  4.3  6 
the  larger  irregular-shaped  particles  have  reasonable 
sources,  one  being  residue  from  the  binder/ inhibitor,  and 
the  second  being  residue  from  the  lab  grade  acetone  used  to 
clean  the  samples. 

D.   COMPARATIVE  EVALUATION  OF  THE  MALVERN  2  600  AND  THE 
LOCALLY  DESIGNED  LASER  DIFFRACTION  SYSTEM 

From  the  outset  of  this  investigation,  a  comparative 

analysis  of  the  two  diffraction  instruments  was  a  primary 

goal.    Some  of  the  limitations  of  the  NPS  system  were 
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obvious,  based  on  a  review  of  the  system  hardware  and 
software. 

The  differences  in  software  and  in  data  reduction 
algorithms  were  significant.  The  NPS  algorithm  is  based  on 
a  ratioing  technique  developed  by  Buchele  [Ref.  10].  This 
method  yields  the  size  parameter  D32  by  the  process  of 
fitting  a  theoretical  curve  for  monomodal  distributions  to  a 
measured  intensity  profile.  The  technique  is  discussed  in 
detail  in  Chapters  II  and  III.  The  MALVERN  Master-Sizer 
operating  system,  by  comparison,  inverts  the  scattered  data 
to  obtain  a  size  distribution.  The  NPS  system  could  be 
upgraded  in  this  area  by  adopting  one  of  the  quadrature 
methods  described  by  Koo  [Ref.  16] . 

The  differences  in  hardware  are  highlighted  by  the 
design  of  the  detector  system  photodiode  arrays.  The 
MALVERN  systems  annular  array  design  allows  for  the 
detection  of  much  weaker  intensities.  The  field  of  view  of 
the  NPS  system  was  also  limited  by  the  array  dimensions. 

A  comparison  was  made  using  the  NPS  laser  system  to 
measure  particle  size  at  the  nozzle  entrance.  Figure  4.37 
is  a  plot  of  comparable  through-motor  measurements  made  by 
the  MALVERN  system  and  the  NPS  system  in  the  motor  section 
just  prior  to  the  nozzle  entrance.  As  can  be  seen  from  the 
plot,  the  NPS  system  underestimated  the  mean  particle  sizes 
at  low  pressures  where  the  distributions  were  multimodal. 
As  discussed  above,  the  NPS  data  reduction  algorithm  assumes 
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a  monomodal  distribution.   The  larger  particles  in  the  bimodal 
distribution  (60-70  microns),  scattered  light  at  angles  below 
the  range  of  the  NPS  detector.   To  further  verify  this  effect, 
the  MALVERN  data  for  the  run  at  161  psig  was  reduced,  after 
inhibiting  data  that  were  outside  the  field  of  view  of  the  NPS 
system.   Figure  4.38  depicts  the  resulting  distribution.   What 
occurred  was  that  an  artificial  mode  appeared  at  24.7  microns. 
Figure  4.39  is  a  measurement  made  at  a  similar  pressure  (175 
psig)  with  the  NPS  system.   A  SMD  of  24  microns  was  obtained. 
Thus,  when  distributions  are  multimodal,  accurate  values  of 
D?_  are  not  obtained  with  present  data  reduction  algorithm 
used  with  the  NPS  apparatus.   Figures  4.40,  4.41  and  4.42  are 
included  for  completeness.   These  figures  represent  the  data 
used  to  construct  Figure  4.37.   At  higher  chamber  pressures 
(400  psig  class) ,  where  the  distributions  became   more  mono- 
modal  and  the  dominant  scattered  light  was  within  the  field 
of  view  of  the  NPS  system,  better  agreement  of  data  existed 
(Figure  4.37)  . 

The  conclusion  is  that  the  NPS  system  functions  well  when 
the  distribution  of  particulates  is  monomodal  and  the  linear 
diode  array  is  placed  such  that  the  majority  of  the  scattered 
light  is  within  the  field  of  view  of  the  detector. 
E.   SOURCES  OF  MEASUREMENT  ERROR  (UNCERTAINTY) 

In  the  course   of   this   investigation   a   number   of 
measurement   errors   were   encountered   from   a  variety   of 
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sources.  Measurements  were  made  in  both  the  motor  chamber 
and  exhaust  environments,  and  each  environment  presented  a 
slightly  different  set  of  difficulties. 

Measurements  made  in  the  motor  environment  were  subject 
to  the  following  sources  of  uncertainty: 

(a)  scattering  by  the  aperture  (window  frame) 

(b)  limited  field  of  view 

(c)  multiple  scattering  effects  due  to  dense 
concentrations 

(d)  surface  reflections  from  the  windows  in  the  motor 

(e)  diameter  of  aperture  limited  scattering  angle  to  a 
maximum  of  6° 

(f)  cleanliness  of  window  surfaces 

(g)  extraneous  light 

(h)  quenching  effects  due  to  chamber  walls,  and  the 
introduction  of  nitrogen  used  as  a  purge  for  the 
window  surfaces 

(i)  vignetting  of  scattered  light,  i.e.,  sample  not 
completely  within  the  range  of  the  focal  length  of 
the  receiving  lens.  In  the  event  vignetting  occurs, 
particle  size  is  biased  to  the  larger  sizes. 

Scattering  by  the  aperture  is  discussed  in  some  detail 
by  Buchele  [Ref .  10]  .  Although  some  scattering  from  this 
source  does  occur,  if  the  particle  concentration  is 
sufficiently  high  (as  in  the  present  case) ,  particle 
scattering  will  dominate. 

The  limited  field  of  view  of  the  NPS  system  did 
significantly  impact  the  scattering  data  obtained.  In  the 
case  of  the  present  investigation,  the  physical  dimensions 
of  the  photodiode  array  did  not  allow  for  scattered  light  to 


91 


be  detected  at  angles  less  than  .01  radians.   Also  in  the 
current  configuration  of  the  NPS  scattering  device,  the  largest 
particle  size  which  could  reasonably  be  detected  was  40  microns 
The  problem  may  be  overcome  in  large  part  by  increasing  the 
dimension  of  the  array.   A  4096  diode  RETICON  array  has  been 
purchased,  and  should  be  installed  in  the  near  future.   This 
will  provide  a  much  enhanced  field  of  view  for  the  NPS  instru- 
ment . 

Multiple  scattering  effects  tend  to  occur  when  the  separat- 
ion distance  between  particles  becomes  small.   Correction 
methods  have  appeared  in  the  literature  [Refs.  27,28],   These 
effects  were  present  in  the  motor  measurements  and  empirical 
corrections  of  the  data  were  required.   A  reduced  combustion 
chamber  width  is  required  to  reduce  this  problem. 

The  problems  associated  with  windowed  surfaces  are  three- 
fold in  nature.   The  first  is  that  some  reflections  are  present 
The  MALVERN  system  in  the  align  mode  allows  sufficient  adjust- 
ment of  the  detector  in  order  to  eliminate  the  effects  due  to 
reflections.   The  second  problem  was  a  result  of  the  physical 
diameter  of  the  window  which  could  pass  light.   This  restricted 
the  maximum  scattering  angle  to  6  .   If  very  small  particles 
had  been  present,  this  would  have  had  a  biasing  effect  on  the 
received  scattering  data.   The  apparent  size  would  have  been 
larger  than  the  actual  size.  The  third  problem  associated  with 
the  windowed  surface  was  to  provide  sufficient  purge  to  keep  these 
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surfaces  clean.  A  trial  and  error  method  of  adjusting  purge 
flow  rates  was  employed.  If  the  windows  had  become 
contaminated  in  a  previous  run,  the  flow  rate  was  increased 
until  the  surfaces  remained  clean  and  usable  data  were 
obtained. 

The  introduction  of  the  purge  nitrogen,  along  with  the 
walls  of  the  motor  chamber,  had  some  undetermined  quenching 
effect  on  the  gas  and  particle  temperatures.  Although  the 
percent  of  gas  flow  comprised  by  the  nitrogen  purge  was 
estimated  to  be  between  3  and  9  percent  for  the  given  range 
of  pressures,  it  was  not  possible  to  determine  what 
influence  the  quenching  effect  had  on  particle  size,  or  the 
extent  to  which  it  actually  occurred. 

Extraneous  light  can  also  be  a  factor.  Buchele  [Ref. 
10]  details  these  effects.  The  ratioing  method  adopted  in 
the  data  reduction  algorithm  for  the  NPS  system  has  the 
desirable  feature  of  eliminating  this  effect.  Reference  10 
points  out  the  fact  that  extraneous  light  is  most  likely  to 
be  a  problem  at  small  angles. 

Vignetting  effects  occur  when  scattered  light  is  not 
captured  by  the  receiving  lens.  This  tends  to  bias  particle 
size  toward  the  larger  size  bands.  The  way  to  overcome  this 
difficulty  is  to  ensure  that  the  field  of  particulates  is 
well  within  the  focal  length  of  the  receiving  lens,  so  that 
all  of  the  light  scattered  by  the  sample  will  be  captured  by 
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the  (Fourier  transform)  lens  and  focused  onto  the  detector 
array. 

The  measurement  errors  that  were  found  to  be  associated 
with  obtaining  scattering  data  in  the  exhaust  are  listed 
below: 

(a)  unstable  refractive  index  of  the  suspending  medium 
which  result  from  thermal  or  density  gradients 

(b)  turbulence  along  the  jet  boundary 

(c)  shocks  associated  with  a  jet  which  is  not  ideally 
expanded 

(d)  the  location  at  which  the  measurements  are  taken 
aft  of  the  nozzle  can  affect  the  measured 
distribution 

(e)  vignetting,  similar  to  that  described  previously 

(f)  extraneous  light  (also  discussed  previously  in 
conjunction  with  motor  measurement  errors) . 

Thermal  gradients  play  a  significant  role  in  measuring 

particle  size  in  hot  gas  jets.   Close  to  the  exit  plane  of 

the  exhaust,  the  jet  was  slightly  smaller  than  the  MALVERN 

laser  beam  diameter.   As  discussed  above,  scattering  data 

were  influenced  by  the  thermal  gradients  along  the  jet 

boundary.   Also  as  discussed  earlier,  the  "KILLDATA"  command 

can  be  used  to  inhibit  data  from  the  inner  rings  where  the 

steering  effect  is  most  prominent.   Buchele  [Ref.  10]  points 

out  that  turbulence  compounds  these   effects.     It  was 

impossible  to  separate  the  effects  of  turbulence  and  thermal 

gradients  on  the  beam  steering.   In  either  case,  so  long  as 

the  user  is  reasonably  assured  that  no  scattering  data  are 

present  as  a  result  of  large  particles  in  the  exhaust  flow, 
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the  "KILLDATA"  command  may  be  used  to  inhibit  data  from  the 
light  steered  to  low  angles. 

The  presence  of  shocks  in  the  exhaust  jet  which  is  not 
ideally  expanded,  is  to  be  expected.  The  density  gradients 
in  the  region  of  a  shock  can  also  adversely  affect  scattered 
data.  If  shocks  are  identified  as  a  problem,  a  schlieren 
may  be  used  to  establish  shock  location  so  that  on 
subseguent  runs  the  apparatus  can  be  positioned  upstream  of 
the  shock. 

The  location  at  which  diffraction  measurements  are  made 
in  a  supersonic  exhaust  can  influence  the  resulting 
distribution  of  particulates  which  are  observed.  The 
further  back  in  the  jet,  the  less  concentrated  the  particle 
field  becomes.  This  may  have  biased  the  measurements  to 
the  smaller  particles  which  were  present. 

The  effects  of  extraneous  light  and  vignetting  have  been 
previously  discussed  and  will  not  be  discussed  further, 
except  to  say  that  the  measurement  and  subtracting  of 
background  light  plays  a  role  in  reducing  the  effect  of 
extraneous  light. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  investigation  have  demonstrated  that 
laser  diffraction  methods  can  be  an  effective  means  for 
determining  particle  sizes  in  solid-propellant  rocket 
motors.  Measurements  were  successfully  made  in  both  the 
motor  chamber  and  exhaust  environments.  Each  environment 
presented  a  different  set  of  difficulties. 

The  MALVERN  and  the  locally  designed  diffraction  system 
produced  data  which  correlated  quite  well  when  the  particle 
distributions  were  monomodal,  and  the  dominant  scattered 
light  was  within  the  field  of  view  of  the  NPS  linear  diode 
array.  When  the  distributions  were  multimodal,  the  NPS 
system  did  not  accurately  measure  the  Sauter  mean  diameter. 

This  problem  was  the  result  of  both  a  limited  field  of 
view  and  the  data  reduction  algorithm.  The  latter  was  based 
upon  monomodal  distributions.  The  system  could  be 
significantly  improved  by  using  larger  arrays  and  by  using 
data  reduction  techniques  which  are  based  on  quadrature 
methods  [Ref.  16]. 

The  NPS  system  will  always  be  limited  to  higher  number 
density  flows  since  the  linear  array  only  intercepts  a  very 
small  fraction  of  the  light  scattered  at  a  given  angle. 

The  MALVERN  2600  has  proven  to  be  well  suited  for  use  in 
solid  propellant  combustion  experiments.   The  beam  diameter, 
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however,  is  larger  than  optimum  for  the  current  application. 
The  acquisition  of  a  new  beam  expander  to  reduce  the  diameter 
of  the  collimated  beam  would  enhance  both  in-situ  motor  meas- 
urements and  exhaust  measurements. 

A  transition  from  the  axisymmetric  motor  to  a  two-dimensional 
motor  would  improve  through-motor  measurements  in  the  following 
ways  : 

(a)  advantage  of  reduced  optical  depth 

(b)  reduction  or  elimination  of  multiple  scattering  effects 
that  may  be  present  in  the  motor  chamber  environment 
due  to  dense  particle  concentrations 

(c)  increased  range  of  scattering  angles 

(d)  allows  operation  to  higher  pressures. 

The  motor  data  suggested  that  ZrC  reacts  to  form  ZrO„  on 
or  near  the  propellant  surface.   Surface  agglomerations  occur- 
red and  were  observed  to  decrease  in  size  as  chamber  pressure 
increased. 

At  the  low  Mach  numbers  in  the  present  motor  the  Zr02 
particles  appeared  to  remain  practically  unchanged  as  they 
passed  down  the  length  of  the  motor  port  at  very  low  pressures. 
At  higher  pressures  (above  250  psig) ,  particle  growth/agglom- 
eration occurred. 

At  low  pressures  the  ZrO-  particles  were  multimodal.   At 
higher  pressures  (  greater  than  350  psi)  the  distributions 
approached  monomodal. 


97 


ZrO~  particles  broke  up  in  the  nozzle  converging  section 
to  a  monomodal  distribution  with  a  D^„  between  15  and  19  microns, 
independent  of  the  particle  size  entering  the  nozzle.   Limited 
data  indicated  that  no  appreciable  change  in  particle  size 
occurred  in  the  supersonic  nozzle  flow.   Aft  of  the  nozzle, 
in  the  expanding  plume,  the  particle  size  was  significantly 
smaller  (2-8  microns).   This  could  be  due  to  particle  breakup, 
or  to  measurement  location  biasing. 

SEM  analysis  of  particles  collected  form  the  converging 
nozzle  wall  showed  particles  with  sizes  in  good  agreement 
with  the  optically  obtained  data. 
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